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ABSTRACT 

We have carried out high-resolution (~3') HI and CO line observations along one- 
dimensional cuts through the Galactic worm GW46.4+5.5. By comparing the HI data 
with IRAS data, we have derived the distributions of iioo excess and ri o excess, which 
are respectively the 100 /im intensity and 100 /zm optical depth in excess of what would 
be expected from HI emission. In two observed regions, we were able to make a de- 
tailed comparison of the infrared excess and the CO emission. We have found that 
r 100 excess has a very good correlation with the integrated intensity of CO emission, 
WcO) but J 100 excess does not. There are two reasons for the poor correlation be- 
tween J 10 o excess and Wqo'- firstly, there are regions with enhanced infrared emissiv- 
ity without CO, and secondly, dust grains associated with molecular gas have a low 
infrared emissivity. In one region, these two factors completely hide the presence of 
molecular gas in the infrared. In the second region, we could identify the area with 
molecular gas, but iioo excess significantly underestimates the column density of molec- 
ular hydrogen because of the second factor mentioned above. We therefore conclude 
that r 10 o excess, rather than 7 10 o excess, is an accurate indicator of molecular content 
along the line of sight. We derive r 100 /iV(H) = (1.00±0.02)x 10~ 5 (10 20 cm" 2 )" 1 , and 
X =iV(H 2 )/Wco— 0.7xl0 20 cur 2 (K km s -1 ) -1 . Our results suggest that Ji o excess 
could still be used to estimate the molecular content if the result is multiplied by a correc- 
tion factor £ c = < I W o/N(H) > m /< Iioo/N(H) > Ha (~ 2 in the second region), which 
accounts for the different infrared emissivities of atomic and molecular gas. We also 
discuss some limitations of this work, which could stem from using single-temperature 
model and the IRAS 60 /im intensity in estimating the dust optical depth along the line 
of sight. 

Subject headings: ISM: individual (GW46.4+5.5) — infrared: ISM: continuum — ISM: 
molecules — radio lines: ISM 
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1. Introduction 

In many studies of the interstellar medium, it is es- 
sential to determine accurately the amount of molec- 
ular gas along the line of sight. However, the most 
abundant molecule, H2, not only has a large excita- 
tion temperature (T ex >509 K), which is not usually 
attainable in a cold interstellar medium, but also has 
very small rotational transition probabilities. In or- 
der to infer the H2 column density iV(H2), therefore, 
we generally rely on indirect methods. The most well- 
known indirect method for estimating the H2 column 
density is to observe other molecules such as CO, 
CS, or NH3. However, the fractional abundance of 
a molecule relative to H2 is unlikely to be uniform 
from one cloud to another, or even sometimes within 
a cloud. In addition, there may be some differences 
among the distributions of individual molecules be- 
cause of their different photochemical properties. In 
the case of CO, which has been widely used as a tracer 
of H2, the abundance varies from cloud to cloud by up 
to three orders of magnitude depending on the visual 
extinction and astrochcmical properties (Scoville & 
Sanders 1987; van Dishoeck & Black 1988; Magnani 
& Onello 1995). 

The infrared all-sky maps produced by the In- 
fraRed Astronomical Satellite (IRAS) mission pre- 
sented a new opportunity to study the distribution 
of interstellar molecular gas. The Galactic radi- 
ation in the far-infrared (far-IR) appears to arise 
mostly from dust grains well-mixed with interstel- 
lar gas (Mathis, Mezger, & Panagia 1983). The 
IRAS 100 (im emission intensity, J100, has been 
found to be tightly correlated with the total col- 
umn density of interstellar gas at \b\ >5°, and the 
100 ^m emissivity per hydrogen nucleus, Iiqq/N{B), 
seems to be fairly uniform in the solar neighborhood 
: Jioo/AT(H)~l MJy sr" 1 (10 20 cm" 2 )- 1 (dc Vries et 
al. 1987; Boulanger & Perault 1988; Heiles, Reach, 
& Koo 1988; Deul & Burton 1990, 1992). Therefore, 
/100 excess, i.e., the IRAS 100 /jm emission in ex- 
cess of what is expected from the HI column density, 
could be used as a tracer of H2 in regions where the 
amount of ionized gas is negligible. Based on such an 
idea, several groups have identified "infrared-excess 
(IR-excess)" clouds in order to investigate the distri- 
bution of molecular gas. However CO observational 
studies on the IR-excess clouds showed that there is a 
substantial difference between the distribution of IR- 
excess clouds and that of CO-emitting clouds (Desert, 



Bazell, & Boulanger 1988; Blitz, Bazell, k Desert 
1990; Hcithausen et al. 1993; Reach, Koo, & Heiles 
1994; Meyerdierks & Heithausen 1996; Reach, Wall, 
& Odegard 1998). This discrepancy between the two 
distributions is in a sense expected, because firstly the 
dust abundance and radiation-field strength vary, and 
secondly the CO abundance is low in diffuse molec- 
ular clouds, and finally the infrared emissivity of the 
dust associated with molecular gas is low. 

Another IRAS method for tracing molecular gas is 
to use 100 fim. optical depth, noo- This might be more 
accurate than using 100 /xm intensity, at least in prin- 
ciple, because optical depth is simply proportional to 
the amount of dust along the line of sight if dust prop- 
erties are the same. Indeed, for several high-latitude 
clouds, 100 /mi (or 60 fim) optical depth has been 
found to be correlated with molecular gas column den- 
sity better than the 100 /im intensity alone (Langcr 
et al. 1989; Snell, Heyer, & Schloerb 1989; Jarrett, 
Dickman, & Hcrbst 1989). 'IRAS' clouds similar to 
the IR-excess clouds have been identified using this 
method, too (Wood et al. 1994). These previous stud- 
ies, which were limited to high latitudes (|6|;>10°), 
however, simply calculated dust optical depth after 
subtracting flat background emission from the 60 and 
100 /im intensity maps and showed that it had a very 
good correlation with the integrated intensities of CO 
isotopes. 

In this paper, we make a detailed comparison 
of 100 /im intensity, 100 /im optical depth, and 
the distribution of CO emission along several one- 
dimensional cuts. Our study differs from most previ- 
ous studies in that we use both high-resolution (~3') 
HI and CO data in the analysis (cf. Reach et al. 
1994). The HI data make it possible to derive accu- 
rate correlations between T100 (or /100) and iV(HI), 
which is used to identify 'rioo-excess' (or ii o-excess) 
regions and to estimate the excess column densities 
of H nuclei. This excess H column density is com- 
pared with the distribution of CO integrated inten- 
sity, WcOy to derive the important conversion fac- 
tor X=N (H2) /Wqo ■ F° r comparison, most previous 
studies (e.g., de Vries et al. 1987; Heithausen & Thad- 
deus 1990) that derived the conversion factor from 
the IRAS data have been based on the comparison 
of -Z100 excess and Wco 7 which may significantly un- 
derestimate X (Magnani & Onello 1995; Reach et al. 
1998). In this paper we compare these two approaches 
in detail and discuss their limitations. 

The object that we study in this paper is the 
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galactic worm GW46.4+5.5. Galactic worms are wig- 
gly, vertical structures which look like worms crawl- 
ing away from the Galactic plane in median-filtered 
HI maps (Heiles 1984; Koo, Heiles, & Reach 1992). 
GW46.4+5.5 is a ^8° long, filamentary structure ex- 
tending vertically from the Galactic plane in both far- 
IR (or HI) and radio continuum emission (Figure 1). 
It is possible that GW46.4+5.5 is the wall of a super- 
shell similar to the North Polar Spur, but at a greater 
distance (Kim, Koo & Heiles 1999). We describe the 
HI and CO line observations in Section 2. In Section 
3 we evaluate I wo and noo excesses using the IRAS 
and HI data, and compare their distributions with 
that of Wco- We discuss the nature of Jioo excess 
in GW46.4+5.5 and some potential problems of Two~ 
excess method in Section 4. The main conclusions are 
summarized in Section 5. 

2. Observations 

HI 21 cm line observations were carried out us- 
ing the 305 m telescope at Arecibo Observatory in 
1990 October. The telescope had a HPBW of 3/3 
and a beam efficiency of 0.8 at 1.4 GHz (Reach, Koo, 
& Heiles 1994). We observed both circular polariza- 
tions simultaneously using two 1024 channel correla- 
tors with 5 MHz bandwidth each, so that the velocity 
resolution was 2.03 km s^ 1 after Hanning smooth- 
ing. Each spectrum was obtained by integrating for 1 
minute using frequency switching. We made a total of 
9 one-dimensional cuts through GW46.4+5.5 at a set 
of constant galactic latitudes. The beam separation 
was 3'. The positions of our one-dimensional cuts are 
listed in Table 1. 

CO J=l— line observations were made in 1994 
February using the 4 m telescope (HPBW=2.'5) at 
Nagoya University in Japan. We obtained 8 one- 
dimensional cuts in the same way as in the HI line 
observations. The observed positions together with 
their reference positions are summarized in Table 1. 
An SIS mixer receiver and a 1664 channel Acousto- 
Optical Spectrometer (AOS) with 40 MHz bandwidth 
were used. The velocity resolution was 0.67 km s _1 
after Gaussian smoothing. The system temperature 
varied in the range 480—700 K during the observing 
sessions depending on weather conditions and eleva- 
tion of the source. Absolute-position switching in- 
stead of frequency switching was used in order to pre- 
vent the contamination of spectra by the atmospheric 
CO emission. The reference positions were checked to 



be free of appreciable (Tp<0.1 K) CO emission. The 
velocity was centered at i>lsr=40 km s _1 and the ve- 
locity coverage was from —19 to +91 km s _1 . The on- 
source integration time was 2 minutes and the typical 
rms noise level was 0.2 K per channel after Gaussian 
smoothing. The intensity scale was calibrated with 
respect to the standard source S140 which was as- 
sumed to have T*=20 K (Yang & Fukui 1992). 

3. Infrared Excess and Molecular Gas 
3.1. HI Gas and Infrared Excess 

Figure 2 shows a sample from our HI spectrum 
at (I, 6)=(45.°10, 4.°00) where the strongest CO line 
emission was detected. There are three peaks at posi- 
tive velocities in the HI spectrum and the peaks in the 
velocity range vlsr— 18— 40 km s" 1 are the compo- 
nents associated with GW46.4+5.5 (Kim et al. 1999). 
The components at negative velocities might be from 
the warped Galactic plane outside of the solar cir- 
cle. The physical and dynamical properties of the 
worm will be discussed in a separate paper (Kim ct 
al. 1999). 

We assume that the HI emission is optically thin 
and compute the HI column density along a given line 
of sight, AT(HI), from 

JV(HI) = 1.822 x 10 18 J T h dv cm" 2 , (1) 

where Tb is the brightness temperature in K and v 
is the velocity in km s _1 . Our assumption seems to 
be valid because the peak temperature of the HI line 
is much lower than the typical spin temperature of 
HI gas, T s =125 K, for most sight lines. The integral 
range is from «lsr=— 100 to +150 km s _1 . Figure 
3a is the plot of the 100 ^im intensity, Iwo, against 
TV (HI) in the observed regions. We can see a strong 
correlation between the two physical parameters. Dif- 
ferent beam sizes of the HI and IRAS observations 
could introduce some scattered data points in the fig- 
ure. But the effect must be small because the beam 
sizes are roughly the same (~3'). The figure shows 
only the data from the regions at b >3° where both 
HI and CO line observations were made. The data at 
lower galactic latitudes are not used because of the 
possible contribution from dust associated with ion- 
ized gas. The open circles, squares, and pentagons 
represent the points with detectable CO emission at 
b = 3°, 4°, and 5°, respectively. The area of each 
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symbol is proportional to the integrated CO intensity. 
The crosses represent the points without detectable 
CO. A least-squares fit excluding the points with de- 
tectable CO yields (7 100 /MJy sr" 1 ) = (1.32±0.02) 
(A(HI)/10 20 cm- 2 ) - 3.37±0.7. Our fitting pro- 
cedure involves two steps: (1) All the crosses are 
fitted to a straight line, and (2) the crosses within 
±2cr deviation from the first fit are fitted to a new 
straight line. The negative zero-intercept could be 
from errors in the subtraction of zodiacal light. The 
derived /ioo/A^HI) agrees with the results of previ- 
ous observational studies. For example, Boulangcr 
et al. (1986), de Vries et al (1987), and Heiles 
ct al. (1988) obtained 1.4±0.3, 1.0±0.4, and -1.3 
MJy sr -1 (10 20 cm -2 ) -1 , respectively. 

We now compare the HI column density and 100 /im 
optical depth. Assuming that the far-IR emission is 
optically thin and that the temperature of dust grains 
is constant along the line of sight, the 60/100 /im color 
temperature, Td, can be derived according to the for- 
mula 
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where c\—hc/k-Q—lAA\ and n is the index in the 
emissivity law, Q a bs(A) — A~™. For the graphite- 
silicate dust grain model of Mathis, Rumpl, & Nors- 
dsieck (1977; hereafter MRN dust model), n lies be- 
tween 1 and 2 (Drainc & Lee 1984). We take n=1.5 
in this paper. The 100 /<m optical depth, ti o, is 
derived from ti o~Iioo / B\(Td) , where B\(T) is the 
Planck function. Figure 3b displays the 100 /im op- 
tical depth versus the HI column density. A strong 
correlation also exists between these two quantities. 
A least-squares fit has been performed in a similar 
way to yield (t 100 /10- 5 ) = (1.00±0.02) (A(HI)/10 20 
cm" 2 ) + 0.01±0.6. The estimated t W0 /N(B.) ratio 
is much smaller than the value < tioo/A(H) > = 6.3 
xl0~ 5 (10 20 cm -2 ) -1 obtained by Boulanger et al. 
(1996) from an analysis of the Cosmic Background 
Explorer (CO BE) mission and the Leiden-Dwingcloo 
HI survey data. This difference between the two val- 
ues cannot be entirely due to the calibration differ- 
ence between the IRAS and CO BE observations. In- 
stead, the difference may be attributed to the effects 
of small, transiently heated dust grains (Langer et al. 
1989) and to the insensitivity of the IRAS observa- 
tions to cold (T d <15 K) dust grains (Snell ct al. 1989; 
Jarrett et al. 1989; Wood et al. 1994). 



If the dust-to-gas ratio and the infrared emissivity 
per hydrogen nucleus are uniform over the atomic and 
molecular regions, the infrared excess can be used to 
estimate the H 2 column density. The infrared excess 
may be determined either from 7 10 o or Tioo using the 
following formulae: 
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where the subscript 'c' indicates the quantity cor- 
rected for offset and the angle bracket indicates an 
average ratio. Since we excluded the data points 
with detectable CO, we may use the derived ra- 
tios as the ratios with respect to the column den- 
sity of hydrogen nuclei, i.e., < iioo/A(H) >=1.32 
MJysr" 1 (10 20 cm" 2 )- 1 and < n o/A(H) >=1.00xl0" 
(10 20 cm -2 ) -1 , respectively. 



3.2. Comparison of Infrared Excess and CO 
Line Intensity 

We have detected CO emission in five regions (Ta- 
ble 2). Figure 2 shows the spectrum of the strongest 
CO line and compares it with the corresponding HI 
spectrum. Note that the velocity range of our CO 
observations was from —19 to +91 km s _1 , so that 
we could not detect the CO emission, if any, asso- 
ciated with the HI gas at negative velocities. It is, 
however, very unlikely that there is appreciable CO 
emission in this outer part of the Galaxy. We, in 
fact, made some test observations with enough ve- 
locity coverage, but could not detect any emission 
(Tp<0.2 K). In any case, we believe that we have de- 
tected the emission from most of the CO gas. The CO 
emission from regions C and D is almost certainly as- 
sociated with GW46.4+5.5, while that from regions 
A and B is probably not because the velocities are 
very different. The CO emission from region E is 
also possibly related to the worm although its cen- 
tral velocity is somewhat lower than the velocity of 
the worm. The physical association of molecular gas 
with GW46.4+5.5 will be discussed in a separate pa- 
per (Kim et al. 1999). As we have mentioned in 
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Section 3.1, the regions at 6=1° (A and B) are not 
included in our analysis because of the complication 
due to the presence of ionized gas. Region E is not 
included because it has only a few observed points. 
We thus limit our subsequent analysis to regions C 
and D. 

Fi gures 4a and 4b compare the values of A^(H 2 ) 
derived from equations (3) and (4) with the inte- 
grated CO line intensity in regions C and D, re- 
spectively. The integral ranges are from v LSR =+20 
to +40 km s _1 for region C and from wlsr=+15 
to +30 km s _1 for region D. These are velocity 
ranges where the CO emission is detected. In re- 
gion C, the -/V(H 2 ) T100 distribution matches very well 
that of WcO) but the 7V(H 2 )j 100 distribution does 
not. The N(K2)i laa distribution differs from that 
of Wco significantly at 6~46°— 47°. In region D, 
on the other hand, both A f (H 2 ) Tl00 and Af(H 2 )j 100 
match very well with Wco although their absolute 
scales differ by about a factor of 2. Hence, noo excess 
appears to trace molecular gas better than 7ioo ex- 
cess, which is also obvious from Figures 3a and 
3b since the points with CO are separated more 
clearly from those without CO in the latter. This is 
more clearly shown in Figure 5, where A^(H 2 ) T100 and 
N(¥L 2 )i mo are compared with Wco ^(H 2 ) T100 is ob- 
viously proportional to Wco whereas iV(H 2 )/ 100 has 
much weaker correlation with Wco- We determined 
A(H 2 ) rioo /Wco=(0.70±0.05)xl0 20 cm- 2 (Kkms" 1 )- 1 
by a least-squares fit. The error quoted represents 
only a statistical error in the fit. This conversion 
factor is much smaller than the estimated value for 
molecular clouds in the Galactic plane, X=(1.8— 4.8) 
x 10 20 cm" 2 (K km s" 1 )" 1 (Scoville & Sanders 1987), 
but comparable to that of high latitude clouds, X 
~ 0.5 x 10 20 cm" 2 (K km s" 1 )" 1 (e.g., de Vries et 
al. 1987; Heithauscn & Thaddeus 1990; Reach et al. 
1998). Since there seems to be a systematic differ- 
ence between the correlations in regions C and D, we 
also determined the slopes of the correlations sepa- 
rately. They are (0.40+0.05) and (0.93+0.06) x 10 20 
cm~ 2 (K km s -1 ) -1 , respectively. This difference 
could be due to variations in the dust-to-gas ratio 
(Magnani & Onello 1995) or due to variations of the 
gas-phase carbon abundance (Heithausen & Mebold 
1989). Further study is needed to address this ques- 
tion. 

As we will discuss in more detail in Section 4.2, 
the correlation between Wco and Af(H 2 )/ 100 is not 
apparent partly because of the existence of a local 



heating source in region C. If we limit the analysis 
to region D, there is a weak correlation between the 
two quantities, N{R 2 )i 10 JWco= (0.59+0.05) xlO 20 
cm -2 (K km s -1 ) -1 . Note that the coefficient of pro- 
portionality is about 1/2 that of the IF / co-A r (H 2 ) T100 
relationship, the significance of which will be dis- 
cussed in next section. 

4. Discussions 

4.1. IR excess as a Molecular Tracer 

4-1.1 Jioo excess as a molecular tracer 
We have found that the J 10 o-excess distribution 
does not correlate well with the CO distribution. The 
discrepancy between the two distributions has also 
been noted in several previous studies (Desert et al. 
1988; Blitz et al. 1990; Heithausen et al. 1993; Reach 
et al. 1994; Meyerdierks & Heithausen 1996; Reach et 
al. 1998). In general, the regions with such discrep- 
ancies are divided into two categories: (a) Iioo-excess 
regions without CO emission, and (b) CO-emitting 
regions without Jioo excess. The sources that fall 
into the first category could be regions with enhanced 
100 /im emissivity. The 100 /xm emissivity is gener- 
ally dependent on the dust-to-gas ratio and dust tem- 
perature. It has been found that the dust-to-gas ratio 
changes with position on an angular scale of ~10° and 
by up to a factor of 4 (Burstein & Heiles 1978; Desert 
et al. 1988). Variation in the dust-to-gas ratios on 
smaller scales is expected, too. The steady-state tem- 
perature of dust grains is determined by their optical 
properties and the intensity of radiation field (e.g., 
Drainc & Anderson 1985). The iioo-excess region at 
Z~46° in region C is due to the strong radiation field 
from a nearby O star (see Section 4.2). Alternatively, 
the sources could be diffuse molecular clouds without 
detectable CO (Blitz et al. 1990; Reach et al. 1994; 
Meyerdierks & Heithausen 1996; Reach et al. 1998). 
This is possible because CO may not be shielded by 
HI, H 2 , and dust from photodissociation, and is un- 
likely to be collisionally excited in the low-density re- 
gion. On the other hand, the sources in the second 
category, i.e., molecular clouds without iioo excess, 
could be present because of the low emissivity of dust 
grains associated with molecular gas. If there were no 
embedded heating source in the molecular cloud, the 
dust temperature drops from the surface to the cen- 
ter of the cloud. According to Mathis et al. (1983), 
for example, the dust temperature drops by about 
3 K from the edge to the center for a cloud with 
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A v =5 mag. The 100 ^.m emissivity of dust grains 
is very sensitive to dust temperature and changes by 
about 40% as dust temperature changes by 1 K over 
the range 18< T d <25 K. At Z~46?3 in region C, the 
low emissivity almost completely counterbalances the 
significant fraction of noo excess from dust grains as- 
sociated with molecular gas. 

The low far-IR emissivity of molecular clouds has 
another important implication for using noo excess as 
a molecular tracer. The noo excess significantly un- 
derestimates the Ei column density. In region D, the 
H 2 column density derived from noo excess is smaller 
than that derived from t wo excess by about a fac- 
tor of 2. This is because in equation (3) we used 
< Iwo/N(H) > obtained from HI line observations. 
If there is a molecular cloud along the line of sight, a 
correct expression for the observed 100 /im emission 
would be 



hoo = < ZiooMH) > HI ^V(HI)+2< JiooMH) > H2 N(R 2 ), 

(5) 

where < 7i 00 /iV(H) > HI and < I 100 /N(H) > H2 are, 
respectively, 100 /jm emissivities per hydrogen nu- 
cleus in atomic and molecular regions. Hence, we 
need to apply a correction factor £ c to iV(H 2 )/ 100 in 
equation (3) in order to obtain an accurate H 2 column 
density: 



£c = < JiooMH) > ffl /< WW(H) >h 2 • (6) 

Since < Ii 00 /N(H) > HI is in general greater than 
< I W0 /N(U) > H2 , £ c >1 and, in GW46.4+5.5, £ c ^2. 
Our correction factor £ c is equivalent to < B v (T d ) > HI 
/< B v (T d ) > H2 used by Magnani & Onello (1995) and 
Reach ct al. (1998). From their study of /ioo-excess 
clouds, Reach et al. (1998), for example, found £ c =3.8 
as an average value for the solar neighborhood. For 
the inner Galaxy, if we use the mean temperature 
of < T d > ffl =21.0±1.0 K and < T d > H2 =19.0±1.0 K 
derived by Sodroski et al. (1994) using the Diffuse 
Infrared Background Experiment (DIRBE) 140 and 
240 /im observations, we obtain £ c =2.1. Our result 
for £ c , however, cannot be compared directly with 
these results because our study is based on the IRAS 
60 yum brightness which is largely affected by the emis- 
sion from small grains. Since B v (T d ) is very sensitive 
to Td, it is necessary to estimate the correction factor 
independently in studying specific regions (Lee, Kim, 
& Koo 1999). 



4-1-2 Tioo excess as a molecular tracer 
According to our results, Tioo excess appears to 
be an accurate indicator of the molecular content 
along the line of sight; it not only traces the presence 
of molecular gas correctly but also increases linearly 
with the amount of molecular gas. There are two rea- 
sons for this. First it separates the intensity peaks 
produced by enhanced infrared emissivity. Second it 
corrects for the low infrared emissivity of molecular 
clouds. In the present work, however, there are two 
limitations as discussed below. 

The first limitation stems from using a mean tem- 
perature in deriving noo along the line of sight. The 
expression for .ZV(H 2 ) in equation (4) is, in fact, an 
exact one if dust properties and the dust-to-gas ratio 
remain the same in atomic and molecular regions. In 
practice, however, there is an error because, in de- 
riving Tioo on the right hand side of equation (4), we 
have used an average color temperature along the line 
of sight (defined in equation (2)) even if we need to 
use different color temperatures in deriving the op- 
tical depths of atomic and molecular regions. It is 
straightforward to estimate the error when there are 
only two dust components with different emissivities 
along the line of sight. If there is atomic gas of optical 
depth rioo (HI) at color temperature Thi and molecu- 
lar gas of optical depth tioo(H 2 ) at Th 2 (<Thi), the ra- 
tio of the optical depth t{ o> derived by using a mean 
color temperature, to the true optical depth tioo is 
given by 
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where r = noo (H 2 ) /noo (HI), At = (c 1 /A 100 )(l/T H2 
— 1/Thi), and the other parameters have the same 
meaning as in Section 3. The derived optical depth is 
always less than the true value and, as expected, the 
difference converges to zero as r approaches zero or 
infinity. Figure 6 shows, as solid lines, how t[ 00 /ti o 
varies in the (r, At) plane. The ratio t( 00 /tioo is af- 
fected more strongly by the warmer dust associated 
with atomic gas in the sense that the error in opti- 
cal depth is larger at noo(H 2 )/noo(HI)=r (>1) than 
that at noo(H 2 )/noo(HI)=l/r for a given value of 
At. Similar figures can be found in Langer et al. 
(1989) and Snell ct al. (1989), but Figure 6 is a gen- 
eralized one. The error in optical depth would induce 



an error m 



The error can be determined 
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by the following formula: 

N ^U 00 r{ 00 /r 100 -l/(l + r) 

^(H 2 ) T100 r/(l + r) • 

The dependence of A(H 2 ) T < oo /A(H 2 ) Tloo on (r, At) 
is displayed by dotted lines in Figure 6. The ra- 
tio A(H 2 ) T ^ oo /A(H 2 ) Tloo i s m general smaller than 
r{ 00 /rioo for a given (r, At) pair. The most strik- 
ing difference between the two quantities is that, as r 
approaches zero, the error in A(H 2 ) Tl00 remains con- 
stant for given At whereas that in tioo tends to zero. 
Since the visual extinction is Ay ~ 1.5 mag at the CO 
peak in region D using -/V(H) / Ay = l.9x 10 21 cm~ 21 mag -1 
(Bohlin, Savage, & Drake 1978), Th 2 is expected to 
be <3 K lower than T m (Mathis et al. 1983). There- 
fore, the error in Af(H 2 ) ri00 may be <,20% in region 
D where rioo(H 2 )/rioo(HI) ~l/2 (See^Figure 6). The 
error in the correction factor £ c is directly related to 
that in A(H 2 ) Tl00 because it is actually derived from 
7V(H 2 ) r;oo /7V(H 2 ) /loo : ^/^ c =7V(H 2 ) T{oo /7V(H 2 ) rioo . 

The second limitation derives from using Iqq in 
estimating the dust optical depth. Interstellar dust 
grains are known to be composed of small, transiently 
heated grains and classical, large grains in thermal 
equilibrium (Draine & Anderson 1985; Desert, Boulanger, 
& Puget 1990). The emission at 60 /im comes from 
both small grains and large grains, while the emission 
at 100 ^m originates mainly from large grains. Desert 
et al. (1990) divided dust grains into three com- 
ponents: Polycyclic Aromatic Hydrocarbins (PAHs), 
Very Small Grains (VSGs), and Big Grains (BGs). In 
their model, the contributions of VSGs are ~60% at 
60 fim and ~15% at 100 /im for the solar neighbor- 
hood. Sodroski et al. (1994) found from an analysis 
of the DIRBE and IRAS data that the average val- 
ues of the contribution of small grains at 60 /im are 
~40% and ^60% in the inner and outer Galaxy, re- 
spectively. IRAS studies of nearby molecular clouds 
suggest that the small grains are present largely in the 
halo surrounding the clouds and that their relative 
abundances vary (Boulanger et al. 1990; Laureijs, 
Clark, & Prusti 1991; Bernard, Boulanger, Puget 
1993). Therefore, ti o that we have derived using 
^60 / Iioo does not represent the usual "optical depth" 
of dust grains. Instead it might be largely affected 
by small grains with varying abundances. However, 
the fact that we do observe a very good correlation 
between noo excess and Wco (or noo and AT (HI)) 
seems to indicate that Tioo is still proportional to the 
amount of dust along the line of sight. Compared 



to iioo, ^ioo carL be interpreted as a quantity cor- 
rected for variation in the abundance of small grains 
as well as variation in the temperature of large grains. 
If one uses far-IR data at wavelengths longer than 
100 Mm, such as the DIRBE data for the 100, 140, 
and 240 fim wavebands, one can resolve this prob- 
lem. In the present days, however, the IRAS data 
have the highest angular resolution among the avail- 
able infrared survey data. 

In the past decade, several groups have found for 
high-latitude clouds that Tioo ( or T 6o) is better corre- 
lated with the gas column density than iioo (Langer 
et al. 1989; Snell et al. 1989; Jarrett et al. 1989; 
Wood et al. 1994). At high latitudes, however, it is 
unlikely that the Tioo-excess distribution is markedly 
different from the iioo-excess distribution, because 
most high-latitude clouds are diffuse and quiescent 
compared to the molecular clouds near the Galactic 
plane. In addition, the low signal-to-noise ratio of 
^60 could introduce a large error in noo- At low lat- 
itudes, on the contrary, the Tioo-excess distribution 
is expected to represent the distribution of molecular 
gas much better than the iioo-excess distribution, be- 
cause the number density of luminous stars increases 
and most molecular clouds contain high-extinction re- 
gions and/or embedded luminous stars. 

4.2. Local Heating Source in GW46.4+5.5 

In the case of region C, the Jioo-excess distribu- 
tion is largely affected by the variation of dust tem- 
perature (Figure 4a). The iioo-excess peak at 6=46? 
appears not to be due to dust grains associated with 
molecular gas but to an increase in dust tempera- 
ture. We found a massive star, LSII +12?3, near 
the position of interest, (I, b) = (46?0, 3?0) (Drilling 
1975). The star was found to be a luminous blue gi- 
ant star, 09 III (Vijapurkar & Drilling 1993). The 
coordinates of the star are (a, /3)i95o=(19:03:14.3, 
12:46:21) or (I, fo)=(45?97, 2?75). An infrared point 
source (IRAS19031+1247) and a radio point source 
with a diffuse extended envelope (Fiirst et al. 1990) 
lie in the vicinity of the star. IRAS19031+1247 seems 
to be an OH/IR star on the basis of its position in 
the (60—25) /im versus (25—12) /im color-color dia- 
gram (cf. Lewis 1994). In order to reveal the physical 
relationship between the two sources and the star, 
further studies are required. We estimate the dis- 
tance of the massive star to be c?=2.2±0.8 kpc, us- 
ing the spectroscopic data mentioned above and the 
UBV photometric data of Drilling (1975), V=10 m 72, 
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B - V=l m OO. Here we adopt M y =-5™3±0 m 7 for 
09 III (Conti ct al. 1983) and R = A v /E B -v=3-3 
because R for the region is likely to be somewhat 
higher than the average value for the Galactic plane 
(Turner 1976). The estimated value is barely in agree- 
ment with the kinematic distance of HI gas associated 
with GW46.4+5.5, d kin =lA kpc (Kim et al. 1999). 

Figure 7 displays the large-scale distribution of the 
ho/hoo ratio around region C. The Ieo/hoo ratio is 
enhanced around the massive star, which is located 
where two arrows meet, and decreases monotonically 
from the star to the position of interest. Assuming 
that Qabs(A) varies as A~ n , the total excess infrared 
luminosity from the I 60 1 1 wo enhanced region can be 
derived from the following formula: 



r CX 

^IR 



= Aad 



j 



< T > Td T d - < T >Td,B T d B 



(9) 

where a is the Stephen-Boltzmann constant, d is 
the distance to the massive star, and the subscript 
B indicates quantities in the absence of the star. 
The Planck-averaged optical depth is defined as < 
r >t= / T X B x (T)d\/ J B\(T)d\ (oc Td 1 - 5 for n=1.5). 
Using d=2.2 kpc and T diB =25.2 K (7 60 //ioo=0.22 
with n—1.5), the total excess infrared luminosity is 
Ljp=6500 L Q over a solid angle of £! s =2.4xlO~ 4 sr. 
This value is an order of magnitude smaller than the 
stellar luminosity, L*=2.2xl0 5 L Q (Panagia 1973). 
Hence, the star has enough luminosity to produce the 
observed infrared luminosity, and the Jioo-excess peak 
at 6=46? in region C is probably due to heating by 
an O-type giant star. 

5. Conclusions 

Infrared (IR) excess, defined as the emission in ex- 
cess of what would be expected from HI emission, has 
been used for estimating the amount of molecular gas 
along the line of sight. However, at the same time, 
there have been molecular, particularly CO, line stud- 
ies showing that the method does not necessarily yield 
reliable results. That is, there are both IR-excess re- 
gions without detectable CO and CO-emitting regions 
without IR excess. Region C in our study is a good 
example. The sources that fall into the first category 
could be either regions with enhanced IR emissivity or 
diffuse molecular clouds without CO emission. In re- 
gion C, it is excess heating from a nearby O star that 
produces an IR-excess region. On the other hand, 



the sources that fall into the second category could 
be molecular clouds with low IR emissivity In region 
C, the low emissivity completely hides the presence 
of a molecular cloud in the infrared. 

Our result shows that we can avoid the above con- 
fusion problem by using the optical-depth excess in- 
stead of IR excess, i.e., the optical depth in excess 
of what would be expected from HI emission. We 
have found a very good correlation between 100 /xm 
optical-depth (tioo) excess and integrated intensity 
of CO emission, Wqo- We derive the conversion fac- 
tor between iV(H 2 ) and W C o, A~0.7xl0 20 for the 
galactic worm GW46.4+5.5. According to our result, 
however, the conversion factors in regions C and D 
differ by about a factor of 2, although both regions 
are parts of the worm. A more complete observation 
is necessary for the study of the variation of the con- 
version factor. 

Another merit of using optical-depth excess, which 
may be more important, is that it gives an accu- 
rate value for the amount of molecular gas. IR ex- 
cess would significantly underestimate iV(H 2 ) if one 
uses an IR emissivity estimated from HI emission, 
< 7ioo/-W(H) > ffl . We have introduced a correc- 
tion factor £ c = < JiooMH) > m /< I W0 /N(R) > H2 , 
which should be applied in order to account for the 
different IR emissivity of molecular clouds. In gen- 
eral £ c > 1 because dust grains in molecular clouds 
are usually colder than dust grains in atomic gas. In 
region D in our study, £ c — 2. We should therefore be 
rather cautious in converting IR excess to iV(H2). 

Two potential problems in using optical-depth ex- 
cess, however, would be the low surface brightness and 
the contribution of small, transiently heated grains. If 
the surface brightness is low, the noise and/or the sys- 
tematic errors produced in the subtraction of zodiacal 
emission may make it impossible to derive an accu- 
rate optical depth. It would be worthwhile to check 
the applicability of the method at high galactic lati- 
tudes. The second problem may be resolved by using 
far-IR data at wavelengths longer than 100 /um. 
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Table 1 

POSITIONS OF ONE-DIMENSIONAL CUTS AND OFF POSITIONS 









/ 


-range 




OFFs a 


Scan Number 


b 


HI 21 cm line 


CO J=l 


-0 line 


(l,b) 


Bl 


1?0 


42?5- 


-48?0 


43?0- 


-45?5 


(45?0, 1?5), (43?5, 2?5) 


B2 


2.0 


43.0- 


-45.7 


43.0- 


-46.0 


(45.0, 1.5), (43.5, 2.5) 


B3 


3.0 


42.5- 


-48.0 


44.0- 


-47.0 


(43.5, 2.5), (45.8, 4.2) 


B4 


4.0 


43.0- 


-47.5 


44.0- 


-46.5 


(47.5, 3.5), (45.8, 4.2) 


B5 


5.0 


42.5- 


-50.0 


45.0- 


-48.0 


(46.5, 5.3), (49.0, 5.5) 


B6 


6.0 


42.5- 


-51.5 


46.0- 


-48.5 


(46.5, 5.3), (49.0, 5.5) 


B7 


7.0 


42.5- 


-51.5 


46.5- 


-47.5 


(45.5, 6.5) 


B8 


8.0 


42.5- 


-51.0 


46.5- 


-47.5 


(48.0, 8.0) 


B9 


9.0 


42.5- 


-51.5 









a OFF positions for CO 3=1— line observations 



Table 2 

REGIONS WITH CO J=l-0 LINE EMISSION 









VLSR a 


J R 


AwpWHM a 


Region 


b 


/-range 


(km s _1 ) 


(K) 


(km s _1 ) 


A 


1?0 


43?50-43?70 


+60 


2.5 


2.1 


B 


1.0 


44.50-44.75 


-10 


2.1 


2.9 


C 


3.0 


45.30-46.65 


+27 


3.0 


7.8 


D 


4.0 


44.20-45.50 


+22 


3.3 


5.6 


E 


5.0 


47.65-47.95 


+ 14 


1.8 


2.5 



a Line parameters at peak position 
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GALACTIC LONGITUDE 



Fig. 1. — 100 /im (grey scale) and 11 cm continuum (contour) images of GW46.4+5.5. The purpose of this figure 
is to show the morphology of GW46.4+5.5 and to show how the observed regions are related to it. We, therefore, 
have removed a smooth background emission from the 100 /im and radio continuum maps using a median filter 
of 6°x2°. GW46.4+5.5 is the far-IR structure extends straight up to 6=7° from the Galactic plane. The regions 
where HI and CO line data were obtained are shown as light and dark solid lines, respectively. For clarity, the 
dark solid lines are shifted in 6-direction from the light ones. The regions with CO line emission are also labeled as 
A— E. Grey scale intensity range is 1—30 MJy sr _1 . Contour levels are 0.03, 0.1, 0.25, 1.0, and 3 K in brightness 
temperature. 
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Fig. 2. — HI 21 cm and CO J=l— line profiles at (I, 6)=(45?10, 4?00) where the strongest CO line emission was 
detected. HI and CO lines were obtained using the Arecibo and Nagoya 4 m telescopes, respectively. The peaks 
between ulsr— 18 and 40 km s -1 in the HI line profile are the components associated with GW46.4+5.5. The peak 
at fLSR — 20 km s _1 has associated molecular gas. 
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Fig. 3. — (a) Jioo vs. AT(HI) and (b) rioo vs. iV(HI). Open circles, squares, and pentagons represent points with 
detectable CO at 6=3°, 4°, and 5°, respectively. The area of each symbol is proportional to Wco- Least-squares 
fits yield (iioo/MJy sr" 1 ) = 1.32(7V(HI)/10 20 cm- 2 ) - 3.37 and (noo/lO" 5 ) = 1.00(iV(HI)/10 20 cm" 2 ) + 0.01, 
respectively. Results of the fits are shown as solid lines. The dotted lines indicate ±2a deviations from the least- 
squares fit. Note that the points with CO are more clearly separated from the points without CO in (b) than in (a). 
Also note that the areas of symbols become larger as the points become further away from the fit in (b), whereas 
they do not in (a). 
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Fig. 4.— Distributions of 7V(H 2 )/ 100 , N(R 2 )t 100 , W C o, and T d at (a) 6=3° and (b) 6=4°. Note that iV(H 2 )/ 100 
differs considerably from either Wqo or ^(H2)n o at b=3° . The difference is partly due to a local heating source, 
ft is also noticeable that dust temperature drops in the regions where CO line emission is detected. The typical 
statistical error bars are given in each panel. 
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Fig. 5. — Comparison of N(H2) TW o (upper panel) and iV(H2)/ 100 (lower panel) with H^co- Filled and open circles 
represent points at 6=3° and 4°, respectively. 7V(H 2 ) Tl00 , but not 7V(H 2 )/ 100 , has a very good correlation with 
W C o- The straight line is a least-squares fit to the data, (TV(H 2 ) T100 /10 20 cm" 2 ) =0.70(W C o/K km s _1 )-0.31. 



15 




jj I I I I I L 

.1 1 10 100 

r=T 100 (H 2 )/T 100 (HI) 



Fig. 6. — Variations of t{ 00 /tioo (solid line) and A r (H 2 ) r | oo /A^(H 2 ) Tloo (dotted line) in (r, At) plane, r = 
tioo(H 2 )/tioo(HI) and At = (ci/Aioo)(1/7h 2 -1/T H i), where Ci/Ai o=144. The t[ qo /t wo ratio is weighted toward 
the warmer dust associated with HI gas, but AT(H2) r ^ oo /./V(H2) Tloo is n °t- The difference between the two quantities 
is larger at smaller r. The cross indicates the position of region D in (r, At) plane. In order to estimate an upper 
limit to the error in iV(H2)ri o m region D, we assume (7hi, Th 2 )=(25.2 K, 22.2 K). Contour levels are 0.1, 0.2, 
•••,0.9. 
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Fig. 7. — Distribution of Ieo/hoo around region C. A luminous blue giant star (LSII+12?3) is located where two 
arrows meet. Note the high intensity ratio around the star. The region where CO line emission is detected is shown 
as a solid line. Grey scale intensity range is 0.20—0.26. The lowest contour level is 0.20 and contour spacing is 0.01. 
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